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Abstract. We show that the steady state weights of all one dimensional exclusion 
models which are mapped to the Zero Range Process (ZRP) can be written in a matrix 
product form, where the required matrices depend only on the steady state weights of 
ZRP. One infinite dimensional representation of these matrices which works for generic 
systems has also been provided. This is in contrast to the usual matrix product ansatz 
which does not always guarantee a solution for the dynamics dependent algebra that 
the matrices need to satisfy. The formulation helps us study the spatial correlations of 
these exclusion processes which are unreasonably difficult to obtain directly from their 
ZRP correspondence. To illustrate this method we reproduce certain known results, 
and then investigate unexplored correlations in some other model systems. 
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One dimensional driven diffusive systems with stochastic dynamics, motivated on 
physical grounds, are studied to understand macroscopic properties of nonequilibrium 
steady states pQ. These systems set a ground for the study of exotic correlations present 
in generic non-equilibrium systems [2]. Due to lack of generic methods analytical studies 
are limited to only those few systems [HI EJ El El E] for which certain specific techniques 
can be applied. The zero range process (ZRP) is one such exactly solvable model which 
displays rich non-equilibrium behaviour, despite of having a simple factorized steady 
state. Several one dimensional systems which are used to model dynamics of avalanches, 
granular systems, interface growth, polymer dynamics, various transport processes, and 
glasses [9] can be mapped to ZRP. 

In ZRP, particles from one box hop to a neighbouring box with a rate that depends 
on the number of particles at the departure box. The boxes are either vacant or 
can accommodate arbitrary number of indistinguishable particles ignoring hard-core 
interaction. Several lattice models [5j El EEl EH] with hard-core interaction (i.e. exclusion 
processes) can be mapped to ZRP when vacant and occupied sites are identified as boxes 
and particles respectively. These exclusion processes have the same steady state weights 
as that of the corresponding ZRPs. However it is practically impossible to calculate the 
spatial correlations on the lattice as they are extremely complicated functions of ZRP 
variables. 

In this Letter, we introduce an exclusion process on a one dimensional ring with 
the dynamics that is equivalent to a generic ZRP. To calculate the spatial correlations 
we provide a formulation where steady state of the lattice system can be expressed in 
a matrix product form; the particles and holes are represented by two non-commuting 
matrices. A specific representation of these matrices, which depends only on the steady 
state weight of the corresponding ZRP, is provided. The formulation can be generically 
used to study spatial correlations of exclusion processes having a ZRP correspondence, 
which is illustrated with a few examples. 

The exclusion process is defined on a one dimensional periodic lattice where the sites 
labeled by i — 1 . . . L can accommodate at best one hard-core particle. Any arbitrary 
configuration of the system is denoted by {s^ where s$ = 0, 1 corresponds to vacant 
and occupied sites respectively. Particles hop to a neighbouring vacant site, say to its 
right, according to the following dynamics, 



where a block of size n, i.e. an uninterrupted sequence of n Is, is denoted by [n]. Here, 
the total number of particles N and thus particle density p = remains conserved. This 
dynamics does not necessarily indicate the presence of a long range interaction in the 
system. In fact several one dimensional driven systems [5j HI], where particles having 
extra degrees of freedom interact through a short range dynamics, can be described 
effectively by Eq. (CQ). 

This exclusion process is equivalent to a ZRP with N particles distributed in 
M = L — N boxes with hop rate u(n), provided 0s and Is are referred to as boxes 




(1) 
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and particles respectively. 

The steady state weight PzRp{ n i, n 2 ■ • •%) of a configuration {n,} in ZRP has a 
product measure for any generic rate u(n) 

M 

PzRp(n 1 ,n 2 . . .n M ) = Y[f(n k ) (2) 

k=l 

with f{n) = nr=i' u (0~ 1 being the weight of a box having n particles [3]. Equation (j2j) 
can be used to calculate the weight of the corresponding configuration of the exclusion 
process as follows. Since the system is both periodic and translationally invariant, the 
weight P({si}) of a configuration {s,} can always be written as P(0[ni]0[n 2 ] . . . 0[tim]) 
where [n&] is the block after k th on the lattice. From this correspondence it is evident 
that 

M 

P({ Si }) = P(0[m] . . . 0[n M }) = PznAri!, . .n M ) = J[ f{n k ). (3) 

k=l 

To obtain the spatial correlation functions of the exclusion process from the 
mapping to ZRP, one needs to express the site variables Sj in terms of the occupation 
numbers {n&}. 

M k 
k=l 1=1 

Clearly, even the two point correlation functions (siSi + j) for arbitrary j are hardly 
possible to calculate. Alternatively let us rewrite the steady state weight of the exclusion 
process on the lattice as 

P({ Sl }) = Tr[X Sl X S2 ...X SL ], (5) 

where non-commuting matrices X Si keep track of the position indices. This form is 
similar to the steady state weights assumed in Matrix Product Ansatz (MPA) [12J. 
However, X s . here satisfy a unique relation [Eq. ([7]) below] which depends only on the 
single box weight f(n) of the corresponding ZRP whereas MPA-matrices need to fulfill 
different algebraic relations depending on the concerned dynamics. 

For simplicity, we use X Q = E and X\ = D. From Eqs. ([3]) and ([5]) we have 

P({ Si }) = Tr [ED ni ED n2 . . . ED nM ] , (6) 

which can be simplified further by taking E = \a)(/3\, where \a) and are yet to be 
determined. This choice of E, along with Eqs. (j3J) and (jSJ), demands that D must satisfy 
the following relation 

(/3|L>» = /(n). (7) 

This equation does not specify |a), (/3\ and D uniquely. For this formulation to work 
it is sufficient to have one suitable choice of \a), (J3\ and D which satisfy the above 
equation. A finite dimensional representation of these matrices is guaranteed only for 
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certain special f(n), e.g. when f(n) is periodic in n, or when f(n) = c n . For generic 
f(n), however, one can use the following infinite-dimensional representation 

oo oo 

£> = $>>(i + l|; (P\ = ( 1 \> and l«)=S/(*-l)|i)> ( 8 ) 

i=i i=i 

where {|z)} is the standard basis in infinite dimension, i.e., \i)k = Si t k. Both (j3\ and 
.D are independent of the systems concerned; (/3\ is a constant vector and D is a shift 
operator which has the following properties 



oo 



Tr[D m ] = 0; and D m ^ a % \i) = ^ a J+m |z), (9) 

i=l i=l 

where {a^} are arbitrary coefficients. It is only \a) which depends on the specific process, 
and can be constructed directly from the set {f{n)}. Thus, the representation (jSJ) works 
for any exclusion process which can be mapped to ZRP; even when f(n) does not have 
any functional dependence on n. 

The first task is to calculate the partition function. This is done in the grand 
canonical ensemble (GCE) by associating the fugacity z with particles (D) resulting in 
Zl(z) = Tr[C L ] where C = zD + E. The weight of the configuration having no vacant 
site is Tr[D L ] = 0. Thus, Zl(z) is the sum of the weights of all other configurations 
with at least one vacant site. 

L L 

Z L (z) = Y, Tl [(zD) n - 1 EC L - n ] = ^((3\C L - n {zD) n ~ l \a). (10) 

n=l n=l 

The partition function can be calculated explicitly for some simple forms of f(n). For 
arbitrary /(n), however, it is better to use the following generating function 

oo 1 



L=l 



where Zl is the coefficient of 7 in the series expansion of Z(z,j). Instead, Z(z,j) 
can be interpreted as the partition function of the system in variable length ensemble 
(VLE). Together z and 7, here, determine the macroscopic variables 

<">-§§ - < 12 » 

In the thermodynamic limit (L) — > 00, where VLE is expected to be equivalent to GCE, 
one of the variables (say 7) becomes a function of the other one (z). Details of this 
equivalence will be discussed elsewhere. 

Any n-point correlation function can be expressed easily in terms of holes denoted 
by Si = 1 — Si. For example (si) = ^[(7-E)(7C) L ~ 1 ], can be written as 

(s l ) = GW(z, 1 ) = l(P\ T -^\a). (13) 

The density p(z,j) = 1 — G^(z, 7), obtained from above, is same as that calculated 
using Eq. (TP2|) in the thermodynamic limit. Since spatial indices do not appear in (sj) 
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it can as well be obtained directly from the ZRP correspondence. The advantage of this 
matrix formulation will be clear in the calculation of higher order correlation functions 
which carry these indices. 

The 2-point correlation functions Gj = (sjSj +J -) are written as 

Gf (*,7) = ^W&^laMzZ^cl") =l j {l- P W\C^\a). (14) 
All other (n + l)-point correlations can be expressed in terms of G^ and G^s, 

n n 

i=i G i=i 

To evaluate Z(£, 7) we need to know (/3|(X — 7C) _1 and (/ — 7Z-D) _1 |a). The later 
can be obtained from Eq. and it can be shown that, 

Here, = Yl n wn f( n ) * s nothing but the single-box grand partition function of the 

corresponding ZRP with fugacity w. Using Eqs. © and (|T5|) in (TTTj) we have 

The advantage of the use of VLE is clear from the above simple and closed form 
expression of the partition function that we obtain for the lattice system. Equation 
(|T5|) can also be used along with Eqs. (|T6|) and ( TT3I) to obtain the average density in 
VLE, 

P ^ l) = FM+>yzF>(>yz) (17) 

(2) 

Higher order correlation functions can be calculated from G m with arbitrary m for 
which one requires 



(/3|G"» = tf"}{p\—L— \a) = { 7 m }- ^ ( "' : J 



'Z-7G' ' 11 J l- 7 F( 7 z)' 

Here, {x n }g(x) denotes the coefficient of x n in the series expansion of g(x). To proceed 
further F(w) needs to be specified. In the following, we illustrate this formalism with 
some examples. 

Our first example is the restricted exclusion process (RASEP) [8] defined on a 
one dimensional periodic lattice with hardcore particles where a particle is allowed to 
move to a neighbouring vacant site, say to its right, if it is preceded by at least ji 
particles. It has been shown, that RASEP and several variations of it undergo a phase 
transition from an absorbing state to an active phase as the density is increased beyond 
a critical value p c = y^-. Exact correlation functions for RASEP are obtained by using 
a generalization of MPA. Let us reproduce these results using the procedure sketched 
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here, as these models can be mapped to a ZRP [7J with u{n) = 0(ra — //), where 0(x) 
is the Heaviside theta function. 

The simplest case /x = 1 corresponds to the dynamics 110 — > 101, which is also 
a special case of some other models studied earlier [13] in different contexts. In the 
following we consider only this case; generalization to \i > 1 is straightforward. Here, 
the single box weight is f(ri) = Q(n) and correspondingly F(w) = jz^- Thus, in this 
case, | at) = (0, 1, 1 . . .) T . Since (/3\a) = 0, weight of any configuration with two or more 
consecutive 0s vanishes identically, as observed in the steady state of RASEP. 

Using Eq. f|T6|) . the partition function of RASEP is given by 

%( Z >V = n 2~Xn v 

(1 — 72; — 7 z z)(l — 72) 

We have checked that {7 L }Z + z L is the GCE partition function obtained earlier p3] 
using MPA. The extra factor z L is the weight of the configuration {s« = 1} which 
does not have a ZRP equivalence. Clearly Z has to be supplemented by Eq. ( 1121) . 
For any given z, the average system size (L) can be varied by varying 7 and the 
thermodynamic limit (L) — > 00 is reached as 7 approaches 7* = ^ z2 + 4z ~ z , i n this 



limit p(z) = 2— — r = 2 1 + \/i+z I ' an d t ne i n + l)-point correlation function 



00 



( Si s i+1 . . . s i+n ) = J2 ^ L Tr[D n+1 C L - n - 1 ] = p^ = p 



P 



Both the density-fugacity relation and the correlation functions, calculated in VLE, are 
same as those obtained earlier in GCE [8]. 

We close the discussion on RASEP with the comment that an alternative finite 
dimensional representation D = I J , |a) = I I and ({3\ = (1 0) which satisfy 



,0 1, 

Eq. (J7J) is identical to the matrices obtained in [8] using MPA. 

As a second example we study another one dimensional system [TT] where particles 
have extra degrees of freedom, namely spin. Each site i of the periodic lattice is either 
vacant or occupied with one + or — kind of particle; correspondingly n = 0, ±1. These 
particles interact through an Ising-like interaction H = — e/4 TjTj + i which is reflected 
in the dynamics: 

+ ; +040+ ; 0-4-0. (18) 

This model, hereafter referred to as EKLS, is an extension of an exactly solvable system 
studied by Katz-Lebowitz-Spohn [6]. It has been shown [TT] that EKLS can be mapped 
to ZRP approximately, ignoring the spin indices, i.e., 0s considered as boxes and both 
+s and — s are considered as particles. The new site variables are now Sj = rf . ZRP 
hop rate for a system having same numbers of + and — particles is 



(\ t v + e (i^ b (t)\ , M 3i;(2 + 6) + 66 /1 + 6 

u(n) = J n = —3- 1 + ; 6(e) = — -77 - s ; v = J + 1 

v 6 \ n J 4(f + e) V 1 — e 

where J n is the current through a block of size n to the leading order in 1/n. The ZRP 
mapping is known to be exact [JS] for the special case e = 0. 
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'0 1 2 3 4 5 



(2) (2) 

Figure 1. (Color online) p and G 2 as a function of z = — hxG\ '. Results from 
numerical simulations (symbols) of EKLS for e = —8/17 with L = 1000 and a = 1 are 
compared with theoretical calculations (solid line). 



As e is increased beyond a critical value 4/5, where b = 2, a phase separation 
transition occurs in EKLS model for large densities. Extensive numerical studies have 
provided evidence [TT] of such a transition. However, to the best of our knowledge, 
the spatial correlations of these models have not been studied yet. Let us use 
the formalism discussed here to calculate these correlations and compare them with 
those obtained from numerical simulations. Although the spatial correlations can 
be found for any arbitrary 6(e), we choose to work with 6 = 1 corresponding to 
e = —8/17 (antiferro magnetic interaction) where it is possible to find simple closed 
form expressions. Then ignoring the prefactor, which only redefines the fugacity, we have 
u(n) = (n + l)/n. Thus f(n) — l/(n + 1) and correspondingly F(w) = —ln(l — w)/w. 
The partition function takes a simple form 



(1 -jz)(z + \n(l-jz))' 



In the thermodynamic limit (L) — > oo, which corresponds to 7* = z x (l — e z ), the 
density is p = 1 — z/(e z — 1). The two-point correlation functions (JHJ) are then, 

J 7 i r°° 

where (A); = A(A + 1) ... (A + I — 1) is the Pochhammer symbol. Explicitly for j = 1,2 

Gf > < 2)= (l-e-*)(2 + Z ) 

Usually the correlation functions in GCE are compared with those in the canonical 

ensemble (having fixed density p = N/L) by expressing z in terms of p. However, here it 

(2) 

is simpler to use z = — In G\ and express p and other correlation functions in terms of 
G\ . In Fig. 1 these results are compared with the 2-point correlation functions obtained 
from the numerical simulations of EKLS with different densities p. An excellent match 
indicates that the mapping of EKLS to ZRP is a good approximation. 
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In summary, we have studied an exclusion process on a one dimensional ring where 
a particle moves to its rightward vacant neighbour with a rate that depends on the size 
of the block to which the particle belongs. The model has a natural correspondence to 
zero range process but positional ordering of the particles is lost in this mapping. Thus 
the spatial correlations of the exclusion process are unreasonably difficult to calculate 
directly from the steady state weights of ZRP. A method is introduced to obtain the 
same, by rewriting the steady state weights of the exclusion process in matrix product 
form. The required matrices are constructed explicitly using only the single box weight 
of the corresponding ZRP. The method is illustrated using the example of RASEP. Since, 
several one dimensional exclusion processes are known to have ZRP correspondence, this 
formulation can be used to study unexplored spatial correlations in such systems. As 
immediate applications, spatial correlations in EKLS model (discussed here), Tonk's gas 
and RASEP with parallel dynamics (to be discussed elsewhere) have been studied. We 
believe that this method is a useful tool for the study of correlations in one dimensional 
non-equilibrium steady states. 

We are grateful to Deepak Dhar for his useful comments on the manuscript. U.B. 
would like to acknowledge thankfully the financial support of the Council of Scientific 
and Industrial Research, India (Grant No. SPM-07/489(0034)/2007). 
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